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Abstract 

The idea of wave mechanics leads naturaUy to assume the well-known 
relation E = hw m. the specific form H = hW, where H is the classi- 
cal Hamiltonian of a particle and W is the dispersion relation of the 
sought-for wave equation. We derive the expression of H in a curved 
spacetime with an electromagnetic field. Then we derive the Dirac 
equation from factorizing the polynomial dispersion equation corre- 
sponding with H. Conversely, summarizing a recent work, we imple- 
ment the geometrical optics approximation into a canonical form of 
the Dirac Lagrangian. Euler-Lagrange equations are thus obtained for 
the amplitude and phase of the wave function. From them, one is led 
to define a 4-velocity field which obeys exactly the classical equation 
of motion. The complete de Broglie relations are then derived exact 
equations. 
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1 Introduction 



1.1 Context of this work 

The long-standing problem of quantum gravity may mean, of course, that 
we should try to better understand gravity and the quantum. More specif- 
ically, it may mean that we should try to better understand the transition 
between the classical and the quantum, especially in a curved spacetime. 
Quantum effects in the classical gravitational field are indeed being observed 
on neutral particles such as neutrons [HI [191 HSj or atoms [HI [13] , with the 
neutrons being spin ^ particles. This together motivates investigating the 
"classical-quantum correspondence" — the correspondence between a classi- 
cal Hamiltonian and a quantum wave equation — for the Dirac equation in a 
curved spacetime. 

1.2 Foregoing results 

In a past work [T], the classical-quantum correspondence was analyzed gener- 
ally from an exact mathematical correspondence, observed by Whitham [20] , 
between a linear wave operator and its dispersion polynomial, and from the 
de Broglie-Schrodinger idea according to which a classical Hamiltonian de- 
scribes the skeleton of a wave pattern. This analysis led later to deriving the 
Dirac equation from the classical Hamiltonian of a relativistic test particle 
in an electromagnetic field or in a curved spacetime [21 [3] . In the latter case, 
this derivation led to two alternative Dirac equations, in which the Dirac 
wave function is a complex four-vector, with the set of the components of 
the Dirac matrices building a (2 1) tensor [21 [3]. (This transformation be- 
haviour may be designated by the acronym "TRD" : tensor representation of 
the Dirac fields.) In order to see if that makes sense physically, the quantum 
mechanics associated with the Dirac equation was then investigated in detail. 

First, it was found [4] that in a Minkowski spacetime in Cartesian coor- 
dinates, the quantum mechanics of the original Dirac equation is the same 
whether, on a coordinate change, the wave function is transformed as a spinor 
and the Dirac matrices are left invariant (which, as is well known, is the stan- 
dard transformation for this case), or if alternatively the TRD transformation 
mode is used. Moreover, the way in which this equivalence was obtained [1] 
makes it obvious that this equivalence holds also with the third transforma- 
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tion mode which can be considered for the Dirac equation, for which the wave 
function is left invariant and the set of the four Dirac matrices is transformed 
as a four-vector. (This is the transformation mode for the standard Dirac 
equation in a gravitational field P [161 [10].) Then it was found that also in a 
general spacetime, the standard equation and the two alternative equations 
based on TRD [5] behave similarly: e.g. the same hermiticity condition for 
the Hamiltonian holds for these equations [S], and similar non-uniqueness 
problems of the Hamiltonian theory occur [5l [6]. 

1.3 Outline of this work 

In this conference paper, we intend to summarize the main part of the recent 
work [8], and to present a few additional results. Those belong to Section [21 
which extends the former derivation of the Dirac equation from the classical 
Hamiltonian of a relativistic test particle [21 [3] to the situation with an elec- 
tromagnetic field and in a curved spacetime. Then, summarizing the main 
results of Ref. [8], Section [3] will go conversely from the Dirac equation to 
the classical motion through the geometrical optics approximation. 

2 From classical motion to Dirac equation 

2.1 Dispersion equation of a wave equation 

Consider a wave equation which is a linear PDE of the second order: 

F^P= ao(X)* + <(X)9^$ + af (X)9^9,^ = 0, (1) 

where X is the position in the space-time, or more generally in the extended 
configuration manifold V of a system of particles. [V has dimension + 1, 
where is the dimension of the configuration manifold.] To be more precise, 
Eq. ([ID is the local expression of the intrinsic differential operator P (which 
acts on smooth sections ip of some vector bundle E with base V) in a local 
chart X '■ ^ ^ (^^) ^^^^ ^ local frame field (cq) on E, with = the 
column matrix made with the components of ip in the frame field (e^), such 
that tp = "^"^Ca in the domain of (cq). The time coordinate is t = x^/c. 

Let us look for "locally plane- wave" solutions: ^(X) = A exp[i6{X)], 
with, at the point Xq G V considered, d^K^^Xo) = 0, where = d^6 
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are the components of the wave covector. Note that Kq = —uj/c, where 
u = —dtO is the angular frequency, and that k = {Kj) is the spatial wave 
covector. (A Latin index j = 1, A^, while a Greek one yU = 0, A^.) 



This leads [H [2] to the dispersion equation: 

Ux{K) = ao{X) + z a't{X)K^ + i^a'^'' {X)K^K, = 0. (2) 

Substituting ^ d^/ i determines the linear operator P uniquely from the 
function (X, K) i— )■ nx(K), which is polynomial in K at fixed X [201 H]. This 
is also true P [3] in the case that G C" with m > 1 in Eq. ([1]), in 

which case the coefficients common to P and IIx are m x m matrices [2], so 
that nx(K) is then an m x m matrix, too. 



2.2 The classical- quantum correspondence 

For any fixed X G V, consider the dispersion equation (|2]): Ilx(K) = 0, here 
assumed scalar (m = 1). This is a polynomial equation for u = —cKq. 
By following smoothly as a function of X G V a particular root, assumed 
real, of this equation, we get a dispersion relation(s): uj = W{]<.;X). (The 
existence of a such real root depending smoothly on X G V is equivalent to 
the existence of a definite wave mode for the PDE ([1]) [20].) As shown by 
Whitham [20] (see also Ref. [T]), the propagation of k obeys a Hamiltonian 
system: 

dAT,- dW dx^' dW , 

On the other hand, according to the wave mechanics inaugurated by 
de Broglie and Schrodinger, a classical Hamiltonian iJ(p,x, t) = H{p;X) 
describes the skeleton of a wave pattern. Then, the wave equation should give 
a dispersion W with the same Hamiltonian trajectories as H. The simplest 
way to get that is to assume that H and W are proportional, H = hW... 
This leads first to the de Broglie relations in traditional form: E = hw and 
p = hk. Then, substituting A'^ d^/i, it leads to the correspondence 
between a classical Hamiltonian and a wave operator. See Refs. [H [2] for 
details. Thus, setting 

= pj (j = 1, ...,X) and Po = (4) 
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we get the de Broglie relations in a condensed form: 



Pf, = hK^ {fi = 0,...,N). (5) 

2.3 The classical-quantum correspondence needs using 
a connection 

The dispersion polynomial nx(K) and the condition dyK^j_{X) = stay in- 
variant only inside any class of coordinate systems connected by "infinitesimally- 
linear" changes [I], i.e., ones such at the point X{{xq)) = X{{xq)) considered: 

0, fi,u,pE{0,...,N}. (6) 



One example [U E] of a such class is constituted by the locally-geodesic 
coordinate systems at X for a pseudo-Riemannian metric g on V, i. e. , 

^^,,p(x) = o, /i,z/,pe {o,...,iV}. (7) 

Specifying, at each point X G V, a class Cx of coordinate systems valid 

in a neighborhood of X, any two of which exchange by a transition map 
satisfying (|6]), is exactly equivalent to choosing a torsionless connection D 
on the tangent bundle TV [3]. Given any such connection, one substitutes 
dfj, ^ Dn into the wave equation ([T]), into the local plane- wave condition 
which rewrites thus D^K^{Xo) = 0, and into the correspondence from the 
dispersion equation (|2]) to the wave equation. That correspondence writes 
thus !■ D^/i. It applies also to the case where "^{X) G C" with m > 1, 
provided that the dispersion equation ([2]) and the wave equation are first- 
order in fact, i.e. ai^'^ = [3]. 

2.4 Hamiltonian of a particle in a curved spacetime 

In a curved spacetime (V, g) with an electromagnetic field of 4-potential 
Vfj,, the world line of a test particle corresponds with an extremum of the 
generally- covariant action integral 

S = -mcds - ^V^dx'', (8) 
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where 

ds'^ = gxpdx^dx^. (9) 
It follows that the motion derives from the following Lagrangian: 

C (x'^, u'") = -mc^g^^uV - (e/c) V^m''^, m'^ = dx" jdi, (10) 

with ^ an arbitrary parameter along the world line of the particle. The 
canonical momenta derived from this Lagrangian are 

= dLld'u!^ = -mc , \ - (e/c)V;. (11) 
^Jgxpu'^u'p 

The Lagrangian (llUp is an "extended Lagrangian" in the sense of Johns [T^ . 
As shown by Johns [12] (Sect. 11.9), we may associate with an extended 
Lagrangian like C a "traditional Lagrangian" L, by setting 

From the latter, we deduce by Legendre transform a "traditional Hamilto- 
nian" i7'(pj, x*^, a;°). The "traditional momenta" pj coincide with the corre- 
sponding "extended momenta" (for /i = 1,...,A^), the latter ones being 
canonically derived from the extended Lagrangian C by Eq. (fTT]) i (Ref. [12], 
Eq. (11.12)). That is, we have 

p, = P, (j = l,...,iV). (13) 

The traditional Hamiltonian is simply (Ref. [12], Eq. (11.14)) 

i^'(p„x^x°) = -Po(x^w"^). (14) 



At this stage, we can specialize the parameter ^ to be the four-length s, 
Eq. dnj. In that case, the vector with components u'^ is the four- velocity, 
u'^ = = dx^/ds. From ([9]), it verifies g^'^u^Uy = 1, as is well known. 
Hence, with ^ = s the canonical momenta (fTTj) become: 

Pp = -mcu^ - {e/c)Vf,. (15) 

Again because g^^u^u^, = 1, they verify the following energy equation: 

g^" (P^ + -V,) (P, + -k) - m'c' = 0. (16) 
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In the expression fll4p of the traditional Hamiltonian if', the time coordinate 
is arbitrary. Let us get H' as function of the same momenta pj and the 
same space coordinates x'', but with the time coordinate t = x^/c. We do 
that directly in Hamilton's equations for H'. We must set: 

Note that Eqs. ( !T3ll and (fTTll are consistent with the definition (jlj). 

2.5 A variant derivation of the Dirac equation 

The dispersion equation associated with the energy equation (fT6|) by wave 
mechanics, i.e., by the de Broglie relations ([5]), is: 

g^'" (hK^ + -V^^ {hK, + ^k) - m^c^ = 0. (18) 

Applying directly the correspondence D^/i to the dispersion equation 

ffTSjl . leads to a specific form of the curved-spacetime Klein-Gordon equation 
[3]. Instead, one may try a factorization with matrix coefficient 0,7'^, etc.: 

lix (K) = [(?^^ (i^^^ + eV^) {K, + eV,)- m"] 1 

=? (« + ^7^^^m)(/3 + <'^.)- (/i=l = c) (19) 

Identifying the coefficients (with noncommutative algebra), and substi- 
tuting Kfj_ D^/i, leads to the Dirac equation [21 13]: 

(z7^ (D^ + leV^) - m)^ = 0, (20) 

with the anticommutation relation, resulting from the identification in Eq. 

7V + 7V = 2^''''l- (21) 



3 From Dirac equation to classical motion 

3.1 General Dirac Lagrangian in a curved spacetime 

The following Lagrangian (density) [B] generalizes the "Dirac Lagrangian" 
(e.g. [HI [S]) valid for the standard Dirac equation in a curved spacetime: 

/ = \ [^1^{D^^) - i^) 7^^ + 2imm] , (22) 
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where \l/ = is the generahzed Dirac adjoint of \1/ = The field X i— )■ 

A{X) designates the hermitizing matrix, such that A'^ = A, (^7'^)^' = A^y'^ 
pn H]. (Here = (M*)^ is the hermitian conjugate of a matrix M.) In 
a curved spacetime, that matrix becomes indeed generally a field [5]. In the 
usual Dirac Lagrangian, the field A is the constant matrix 7^*^, where (7^") 
is a set of constant "fiat" Dirac matrices, i.e., ones obeying the anticommu- 
tation relation ( 12 ip with the Minkowski metric 77°^ — provided the set (7''°) 
is chosen such that 7^" be a hermitizing matrix for that set. 

The Euler-Lagrange equation for the Lagrangian (l22l) is the generalized 
Dirac equation [5], |6] : 

rD^^ = -zm^ - ^A-\D^{A^^^))^. (23) 

This coincides with the usual form iff D^(A7'') = 0. That is always the case 
[5] for the standard Dirac equation in a curved spacetime (the "Dirac-Fock- 
Weyl" or DFW equation). In Eqs. fl22|) and fl23|) . the covariant derivatives 
Dfj^ correspond to an arbitrary connection D defined on the complex vector 
bundle E, in which the Dirac wave function ^ is living. That vector bundle 
is assigned to be a "spinor bundle", i.e. essentially, one for which it exists a 
global field 7 of Dirac matrices, consistent with the anticommutation relation 
(ED). See Ref. [7j for details. 

3.2 Local similarity (or gauge) transformations 

Admissible coefficient fields (7^^, A) for the general Dirac equation (l23ll are 
ones such that the anticommutation relation ( 12T1) is satisfied and that the 
field A is hermitizing. Given any local similarity transformation S : X 1— >• 
S{X) G GL(4, C), depending smoothly on X G V, other admissible coefficient 
fields are 

Y = S-^YS (;u = 0,...,3), A = S^AS, (24) 

in the sense that the anticommutation relation (!2T|) remains satisfied [in the 
same spacetime (V, gi) !] with the new field of Dirac matrices 7^, and more- 

^ The connection D and the covariant derivatives extend as usual to the dual 
bundle E° of E, and to tensor products such as E (g) E°. Moreover, on a tensor product 
such as TV (g) E (g) E° (which is the relevant bundle for the Dirac matrices [7| ) , the relevant 
connection is got from considering the Levi-Civita connection on the component bundle 
TV [HIT]. 



8 



over the matrix field A is a hermitizing matrix field for 7^ [S]. 

The relevant Hilbert space scalar product is [5] 

I $) = J ^U7°$v^ d^^. (25) 

It transforms isometrically under the gauge transformation (1241) , if one trans- 
forms the wave function according to ^ = S^^"^ [Ej. 

The Dirac equation (123|) is covariant under the similarity (l24l) . if the 
connection matrices F^, such that = + T^, change thus [lOl |6]: 

f, = S-'T,S + S-\d,S). (26) 

For the DFW equation, the gauge transformation ( 1241) is restricted to be- 
long to the spin group: VX S{X) G Spin(l,3). Then, the relation (126|) is 
automatically satisfied [Ej. 

3.3 Reduction of the Dirac equation to canonical form 

If D^{A'y^) = and the F^'s are zero, the Dirac equation (123|1 writes 

= -im^. (27) 

Theorem 1 |j8]. In the neighborhood of any event X , the Dirac equation 
[23\] can he put into the canonical form [2l\ ) by a local similarity transfor- 
mation. 

Outline of the proof: i) By Theorem 3 of Ref. [5] [Sect. 3.4, Eq. (54)], a 
similarity T defined by Eq. [with T instead of 5*] brings the general 

Dirac equation (l23l) to the "normal" form {Dfj^{A'-f^) =0), iff 

A^D^T = -{1/2)[D^{AY)]T. (28) 

ii) By Theorem 2 of Ref. [7] [Sect. 6.2, Eq. (65)], a similarity S defined by 
Eq. (124|) and such that Eq. (l26l) is satisfied, brings a normal Dirac equation 
to the canonical form (!27|) . iff 

AYdf^S = -AYT^S. (29) 
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iii) Due to the hermitizing character of matrix A, and due to the fact that, 
by construction, the Hermitian matrix A'j^ which defines the scalar product 
fl25|) is positive definite [5], both f l28|l and f l29|) are symmetric hyperbohc 
systems. □ 



3.4 Geometrical optics approximation of Dirac Lagrangian 

The Lagrangian for the Dirac equation in an electromagnetic (e.m.) field is 
got by substituting + ||V^, where e is the particle's electric charge. 

It follows that, for the canonical Dirac equation ( 127|1 . the Lagrangian in an 
e.m. field is 



I 



— ihc 



2imc 

h 



(30) 



0, where is the covariant derivative obtained on the 



with V^{A-i^) 

relevant tensor product bundle (see Footnote 1) if the connection matrices 
of the connection D are zero when D is acting specifically on the bundle E 
[8]. To implement the geometrical optics approximation following Whitham 
[20] . we set first 

^ = Xe'\ (31) 

where x = x(^) is also a complex wave function with four components, 
and 9 = 0{X) is a real phase. We assume then that x is slowly changing 
compared to the rapidly changing phase 6. That is, the geometrical optics 
approximation consists in assuming that 

d>.X « {d,9)x. (32) 
Substituting (l3Til into the Lagrangian (l30l) with this approximation, yields 

(33) 



I' = [[-^d^O - -V^j x^A-ff'x - mcx^Ax 



The Euler-Lagrange equations are then [8]: 

(^-hd^e - ^V^^ A-ff'x = mcAx, 
d^. {cV^ X^A^'^x) = 0. 



(34) 
(35) 
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3.5 Classical trajectories 

Theorem 2 jS]. From \1/ = xe*^, define a four-vector field and a scalar 
field J thus: 

= -—d,e - -^V,, (36) 
mc mc' 

u'' = g'"'u,, (37) 
J = c x^Ax. (38) 

Then the Euler-Lagrange equations (34\ ) imply 

V,{Jun = 0, (39) 
= 1, (40) 
V^u^ - "^uU^, = -{e/mc^) F^^. (41) 

The first equation is the conservation of the probability current. The two 
last equations imply the classical equation of motion for a test particle in an 
electromagnetic field in a curved spacetime. 



3.6 De Broglie relations 

The canonical momenta of a classical particle are given by Eq. f|T5|l : 

= -mcu^ - {e/c)V^. (42) 



On the other hand, following a Dirac quantum particle in the geometrical 
approximation, we were led to define a 4-velocity field from the phase 6 
of the wave function, Eq. (1361) : 

= -—d,e - -^V^. (43) 
mc mc^ 

We saw that the field obeys exactly the classical equations of motion, 
which are Hamiltonian equations for which the canonical momenta are given 
by Eq. But, remembering the definition = d^9 of the wave covector 

from the phase 6 of the wave function, we rewrite Eq. (H3l) as 

-mcu^ - ie/c)V^ = KK^. (44) 

That is, we derive exactly the de Broglie relations 
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4 Conclusion 



The Dirac equation in a curved spacetime with electromagnetic field may be 
"derived" from the classical Hamiltonian H of a. relativistic test particle. One 
has to postulate H = hW (i.e., E = hw), where W is the dispersion relation 
of the sought-for wave equation. Then one factorizes the obtained dispersion 
polynomial. 

Conversely, to describe "wave packet" motion, we implemented the geo- 
metrical optics approximation into a canonical form of the Dirac Lagrangian. 
From the equations obtained thus for the amplitude and phase of the wave 
function, one is led to define a 4-velocity field. This obeys exactly the clas- 
sical equation of motion. 

The de Broglie relations = ^A'^ are then derived exact equations. 
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